Break up of Rydberg superatoms via dipole-dipole interactions 
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We investigate resonant dipole-dipole interactions between two "superatoms" of different angular 
momentum, consisting of two Rydberg-blockaded atom clouds where each of them carries initially 
a coherently shared single excitation. We demonstrate that the dipole-dipole interaction breaks up 
the superatoms by removing the excitations from the clouds. The dynamics is akin to an ensemble 
average over systems where only one atom per cloud participates in entangled motion and excitation 
transfer. We also consider a modified setup, in which instead of a single shared Rydberg excitation 
each atom carries a small admixture of an excitation via off-resonant laser dressing. In this case, 
we observe collective motion for the entire cloud. This collective behavior constitutes a further step 
towards the implementation of entangled motional states. Moreover, our findings should facilitate 
the experimental realization of adiabatic exciton transport in Rydberg systems by replacing single 
sites with atom clouds. 
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I. INTRODUCTION 

Rydberg atoms have several remarkable properties, 
such as long lifetimes, large polarizability and strong 
long-range interactions. These properties have turned 
Rydberg atoms into versatile tools for quantum informa- 
tion [IH5], nonlinear quantum optics 6-8 or cavity quan- 
tum electrodynamics [91411]. Recently, also the potential 
of Rydberg atoms for quantum transport in atomic ag- 
gregates has been demonstrated [TMT6] . The underlying 
physical mechanism is closely related to the one of exci- 
tation transport in molecular systems [T3] E3 [TS] . 

So far, only single atoms have been considered as build- 
ing blocks for excitation transport in Rydberg aggre- 
gates. An experimental implementation would be facili- 
tated if the single atom sites could be replaced by atom 
clouds. This has motivated us to explore the possibilities 
for extending the existing Rydberg transport schemes us- 
ing atom clouds as sites. 

When several atoms are brought together in clouds, 
they exhibit an effect known as Rydberg blockade 
[28] . At small interatomic distances, the van-der-Waals 
interactions £/vdw = C 6 /r 6 (where C 6 is a state- 
dependent interaction constant and r the interatomic dis- 
tance) become large. These interactions lead to an en- 
ergy offset of all many-particle states with more than 
one Rydberg excitation with respect to the energy of 
the states with just a single excitation. As a conse- 
quence, although irradiating an entire atomic cloud, a 
laser can create at most one Rydberg excitation reso- 
nantly, within a radius at which this energy offset is 
larger than the laser linewidth. The latter is limited 
by the Rabi frequency f2 of the driven transition, which 
gives an estimate of the minimum interatomic distance 
between two Rydberg excitations, known as the blockade 
radius ru ~ (Cg/M^) 1 / 6 . The interplay of the Rydberg 
blockade with quantum transport has appealing new as- 
pects which we investigate in the following. 

We concentrate on a system of two atom clouds, well 



separated in space. A similar arrangement was stud- 
ied in Ref. [39] in the limit of a frozen Rydberg gas. 
Here we choose the spatial extension a of each cloud 
to be smaller than the blockade radius r^i, while the 
inter-cloud distance L is significantly larger than rbi, see 
Fig. [T] For such parameters, the electron dynamics is re- 
stricted to a single Rydberg excitation per cloud, and we 
consider a situation in which each of the two Rydberg- 
blockaded clouds is initially prepared in a coherent col- 
lective excitation. Such states can nowadays be exper- 
imentally created [24] and are sometimes referred to as 
superatoms |30[ 131] . They are a coherent superposition 
of states where all atoms within the cloud but one are in 
the ground state and one atom is in a Rydberg state 
Here, v and I denote the principal and angular momen- 
tum quantum numbers of the Rydberg state, respectively. 
We shall further imply with the term superatom, that all 
the atoms taking part in the superposition of electronic 
states share almost identical spatial probability distribu- 
tions. 



We will analyze the resonant dipole-dipole interac- 
tions [3"2"r[3"l] between two superatoms of different angular 
momentum, with the setup and the model described in 
section [II] Then, we will discuss the electronic excitation 
transfer in the limit of frozen atomic positions in sec- 
tion |III| which already singles out some interesting fea- 
tures of the two interacting superatoms. Subsequently, 
we will turn to the effects of the electronic excitations on 
the motion of the atoms in section |IV| Finally, in sec- 
tion [V] we will investigate how the character of the dy- 
namics changes, if Rydberg-dressed [35HM| ground state 
atoms are considered instead of an actual Rydberg exci- 
tation. The paper ends with a summary and conclusion 
in section IVI1 
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II. THE SETUP 

Wc denote the two atom clouds as A and B with Na 
and Nb atoms in the respective cloud, while N — Na + 
Nb is the total number of atoms. Each atom is labeled 
uniquely, such that atoms with numbers I.-Na belong to 
cloud A while (Na + l)..N refers to atoms in cloud B, 
summarized in the index sets 

A=[l,..,N A ], B = [N A + 1,..,N}. (1) 

We consider one spatial dimension along the separation 
of the clouds, and take three electronic states into ac- 
count: The ground state \g) and two Rydberg states \vs) 
and \vp) . Since the principal quantum number v of the 
Rydberg states is kept fixed it will be omitted in the fol- 
lowing. With these restricitions, the Hamiltonian for our 
system reads 

N h 2 

H(r 1 ,...,r N ) = -^2^V 2 ri +H e} (r 1 ,... ) r N ), (2) 

i=l 

where M is the mass of a single atom, the position of 
the «-th atom and H c \ the electronic Hamiltonian describ- 
ing dipole-dipole-interactions between the atoms. Gen- 
erally, any of the TV atoms can be in either of the three 
states |<7), \s), \p). However, the Rydberg blockade and 
the binary character of the dipole-dipole interactions sig- 
nificantly reduce the dimensionality of our problem. The 
former excludes states with more than one Rydberg ex- 
citation with the same angular quantum number within 
one cloud, while the latter excludes states in which both 
the s and the p excitation are localized within one cloud. 
We define electronic states 

Kmn) = \ggg~s~ggg~p~ggg) (3) 

where the n-th atom is in the Rydberg state \s), the fri- 
th atom in the Rydberg state \p) and all others in the 
ground state, and impose and additional constraint 

(n € A and to € B) OR (to g A and n <E B) , (4) 

as sketched in Fig. [I] The states of Eq. ([!§ con- 
strained according to Eq. Q thus constitute a basis in 
the 2iV J 4iVB-dimensional electronic subspace. We assume 
that the dipole-dipole interaction of strength Vq 

Vq 

Vnm( r nm) = ~~ 3 (5) 

between the two excited atoms n and to depends only on 
the interatomic distance r nm — \r n — r m \ neglecting possi- 
ble orientation effects. This is a sufficient description for 
certain geometries and selected states [16j [29l |39] . Using 
Eq. ([5]) the matrix elements of the electronic Hamiltonian 
H c \ from Eq. ([2| read 

(^*nm I H c \ |7T n ' m /) — Vnmfinm' ^mn' ■ (6) 




FIG. 1: Color online: Sketch of two interacting Rydberg 
blockaded clouds (not drawn to scale). The width a of 
the Gaussian distribution of the atoms within each cloud is 
smaller than the blockade radius rbi, while the distance L be- 
tween the clouds is larger. The coloring (shading) illustrates 
the electronic basis states j7r nm ) introduced in Eq. Q. 



TABLE I: Matrix elements of the electronic Hamiltonian H E \ 
in the basis |7r nm ) for the example of Na = Nb = 2. 
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As a consequence of the Rydberg blockade, the electronic 
Hamiltonian in the matrix representation, Eq. ([6]), has a 
block-diagonal structure, as can be seen in Table [TJ This 
has drastic consequences for the electron and atom dy- 
namics: (i) The excitation exchange between the clouds 
can be resembled by an average over single atom pairs 
and (ii) the force resulting from the dipole-dipole interac- 
tions only acts on one atom rather than the whole cloud. 
We will discuss these effects in sections |III| and |IV| re- 
spectively. 

The full wave function can be written as an expansion 
in the electronic basis |7r nm ) 

) (7) 

with coefficients <^ nm (R, t) which depend on time t and 
the atomic positions r iy summarized in the vector R = 
(n, .., rjv)- The probability to find an s-excitation on 
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atom n and a p-excitation on atom m is given by 

p nm (t) = J d^Rl^CiM)! 2 . (8) 

Often the initial state of the system can be written as a 
direct product 



Xsp , 



(9) 



where | i/j™ ) is the initial electronic state and | x S p' ) the 
initial wave function of the atoms. The latter is taken 
as a product of one-dimensional Gaussians in position 
space, centered around the respective centers ta, rg of 
the two clouds, 



N A 
n=l 



JV 



n=N A + l 



r B ), (10) 



where 



X£(x) 



(Tr^-^expC-lxl 2 ^ 2 ) (11) 



denote a Gaussian in T> dimensions. The width a is taken 
such that a < r\,i and a <C L, where L = \ta — rs\ is the 
distance between the clouds. This choice of the initial 
spatial wave function models an experimental prepara- 
tion of the atoms in the lowest oscillator states of two har- 
monic traps. Our initial electronic wave function | ) 
is taken to represent a superatom state in both clouds. 
Specifically, all atoms in cloud A coherently share the Ry- 
dberg s-excitation and all atoms in cloud B coherently 
share the Rydberg p-excitation: 



1 

VNa~N~e 



E 



(12) 



neA, meB 



Note that the position-dependent dipole-dipole poten- 
tial entangles the electronic interactions and atomic mo- 
tion [TH - TIBl !"IU] , so that for t > the direct product form 
of the wave function will not persist. 



III. DYNAMICS IN THE ELECTRONIC 
SUBSPACE 

Before tackling the intricate interplay of electronic ex- 
citation and atom motion, we first study purely electronic 
dephasing of population oscillations arising from the ini- 
tial state, Eq. (12). The observable of interest is P|, the 
s-population in cloud A, which can be defined in terms 
of the occupations of the electronic states (cf. Eq. 



m) = 



P-nm (^) ■ 



(13) 



For our initial state (Eq. (12|) we have P%{t = 0) = 1. 
Due to the interactions (Eq. (pH)), the s-population will 



migrate from cloud A to cloud B, while p-population 
will migrate from cloud B to cloud A, resulting in Rabi- 
oscillations of the quantity P\{t). The oscillations de- 
phase due to the width of the initial spatial wave func- 
tions | Xsp ) ■ I n this section, we will show that this be- 
havior can be understood already within a semiclassical 
picture. Moreover, we will also show that the same de- 
phasing can be equivalently obtained as an effect of dis- 
order by statistically avering over fixed point particles. 

For the following semiclassical treatment, it is conve- 
nient to represent the electronic subspace in the basis of 
the eigenstates ^(R)) of the electronic Hamiltonian, 
rather than in the basis |7r nm ) introduced in Eq. (J3j) . 
These states satisfy 

tf cl (R)MR)) = U k {R)\ Vk {R)), (14) 

and the wave function is expanded as 

|tf(R,t))=^ fe (R,t)K(R)). (15) 



We call the eigenenergies £4(R) (adiabatic) potential en- 
ergy surfaces and <^fc(R, t) spatial wave functions in the 
adiabatic representation [T?J HT] . The index k runs over 
the 2NaNb dimensions of the electronic subspace. For 
an electronic Hamiltonian with the structure of Eq. ^ 
(see also Tabic [IJ, the mapping of the basis functions 
from the expansion in Eq. |7| , |7T ram ), onto the adi- 
abatic basis |(^fc(R)), Eq. (15), is particularly simple. 
The eigenvalues U k , k = 1, ..,2NaNb of the electronic 
Hamiltonian can be grouped into NaNb pairs, each of 
which can be written in terms of the Hamiltonian ma- 
trix elements as J7± m (R) = ±V nm (r nm ), cf. Eq. 
Similarly, the corresponding eigenvectors have the form 
|7r m „) ± |7r mn ))/v2. Therefore, we can rewrite 



the adiabatic expansion from Eq. (15) as 



(R,<))= (&m( R >*)K 



nm/ 1 Tnm\ ' /\rnm 



neA.meB 



(16) 

The equations of motion for the coefficients ^nm are then 
obtained from the time-dependent Schrodinger equation 
and read 



ih-4>± m (R,t) = - ^ — V^±„(R,t) 

i—1 



+ UL(r nm )^ m (Rt). (17) 



Since the above equations decouple for different atom 
pairs nm and surfaces ±, we can solve them separately 
and will often suppress the indices nm in the remain- 
der of this section. In the following two subsections, 
we will introduce a semiclassical propagator and show 
that th e latt er yields dephasing oscillations of P% (sub- 
section III A ) , which agree with a full quantum mechani- 
cal calculation. In subsection |III B| we demonstrate that 
exactly the same dephasing oscillations arise even if the 
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atoms are treated as fixed classical point particles. To- 
gether, subsections HI A| and |III B rigorously justify the 
use of the frequently employed so-called frozen gas ap- 
proximation for our scenario. The reader content with 
its applicability may thus directly proceed to the semi- 
classical expression for P%{t) in Eq. (31). 



A. Semiclassical treatment 

Now we introduce the semiclassical approach. We 
choose the Herman-Kluk propagator K HK in the initial 
value representation [42-45], that achieves 



k± (R",t) = I dR'K HK (R",t;B.',t = 0) ± 4> ± (R',t=0) 



The propagation is carried out in terms of multiple classi- 
cal initial value solutions and trajectory averaging. Such 
an approach constitutes a useful numerical method; Here, 
however, we shall use it for some simplifying theoretical 
arguments. 

The propagator reads 

K HK (K",t;Ii',t = 0) ± = 
d N p , d N q' 



R"U)C(q t)P i,q',p' 



x exp ( ^S ± [q t ,Pt 



2h 



(qtPt-qV) )(z'\R'). (19) 



Here, | z' ) are coherent states centered on position q' and 
momentum p', 



X#(R-q')exp 



(20) 



p|/2M on these timescales is negligible compared to the 
potential energy E c p = U ± {q t ). The timescale is essen- 
tially defined by the ratio {L/a) and the inverse uncer- 
tainty of the potential energy, i.e. we restrict ourselves 
to times 

t «*AE p W 

where the uncertainty of the potential energy is estimated 
as 



dr r , 



3V <7 



L 4 



(23) 



With these conditions, as shown in Appendix |XJ the 
Herman-Kluk propagator becomes 



(18) K HK (R",t;R',t = Q) ± = exp 



-§^(R') 



<J(R' - R" 



(24) 

Hence, the time-propagated adiabatic spatial wave func- 
tion is given by 



: (H",t) = exp 



it 



(R",i = 0). (25) 



This result is equivalent to the solution of the Eqs. (17) 



if, in the latter, the kinetic energy operator is ignored. 

Now recall that the adiabatic surfaces C/ ± (R") = 
U^ m (H") not only depend on just two coordinates r n , 
r m , but actually solely on the relative coordinate r nm = 
T n — T m . Further, as already mentioned, the equations 
of motion for the adiabatic wave functions decouple for 
different atom pairs. We can thus without loss of gener- 
ality restrict our treatment to a single atom pair nm und 
transform the wave function onto relative and center of 
mass coordinates: 



R 



CM 
nm 



(26) 



with Xn defined in Eq. ( 11 ). The width <jq of these states 
is just a numerical parameter, which is adapted to the 
given problem. For simplicity, in the following we set it 
to be equal to the width of the initial wave function de- 
fined in Eq. (10), i.e. <7q = a. The expressions q t , p t 



are trajectories evolved with the classical Hamiltonian 
H c = p 2 /2M + U^(qi) from the initial conditions q', p'. 
We denote the classical action along these trajectories by 
S±[q t , Pt ] = J[p 2 /2M-U ± (q t )]dt. The time-dependent 
coherent states centered around the phase space trajec- 
tories qt , pt are labeled by \z t ). Finally, the prefactor 
C is defined by the trajectories as 



C 



dp t <9qt 
dp' dq' 



ih dqt 
a 2 dp' 



<J 2 dp t 
ih dq' 



(21) 



In Appendix [XJ we show step by step that under cer- 
tain conditions, the form of the propagator (19) can be 
substantially simplified. The conditions are (i) L a, 
(ii) the dynamics is considered only on sufficiently short 
timescales, and (iii) the classical kinetic energy E^ — 



After separating the center of mass motion, the dynam- 
ics becomes effectively one-dimensional. From the initial 
state (Eqs. ( |9|To| ), we obtain 

,t = 0) = ^ 7 =xf' r (rn m -L), (27) 



Kim( r nm- • 



with x^ a as defined in Eq. |ll]). The width y/2a results 
from the transformation onto the relative coordinate, 
while the additional \j\[2 factor follows from our nor- 
malization convention in the adiabatic expansion. The 
Herman-Kluk propagation p4|) therefore yields 



4 > nm( r nm, t) 



■- exp 



it 



r 3 



Al V n 



L). (28) 



The s-population in cloud A in the semiclassical picture 
Pg C A for a single atom pair nm can be extracted from the 
spatial wave function in the adiabatic representation as 

1 r 2 



(29) 



5 



Inserting the Herman-Kluk-propagated wave functions 
from Eq. ( 28 1 leads to the expression 



raw 



1/2 



dL ' C0S2 (j^) 



{L'-Lf 



2a 2 
(30) 

This integral cannot be solved analytically exactly, but 
by substituting rj = L — L' and Taylor-expanding [L + 
7]}~ ~ L~ 3 (l — 3rj/L) one finds the following approxi- 
mate solution for the s-population in cloud A: 



raw 



1 + cos(wi) exp 




(31) 



with bj = (2Vq /L 3 )/h. We thus observe superexpo- 
nentially dephasing Rabi-oscillations already within the 
semi-classical picture, where the dephasing time is deter- 
mined by the ratio cr/L: 



7"dcph 



V2L V2hL 4 



3<juj 



6crV 



(32) 



This dephasing time is always within the timescale 
Tdoph <C (L / a) 2 hL 3 / (3Vq) as long as a <C L, which 
is consistent with the conditions initially imposed for the 



semiclassical treatment, cf. Eq. (22 1 



B. Classical point particle limit 

Let us finally remark that it is, in fact, possible, to 
arrive at the very same expression by an even cruder 
approximation, namely by simply describing the atoms 
as fixed point particles and averaging over their positions. 
In such a picture, there is no dynamics in position space 
and the expansion of the wave function in terms of the 
basis |7r„ m ) as given in Eq. can be simply written as 



(33) 



The expansion coefficients f nm evolve according to 
d 

tk-Q^fnrn(t) = V nm (r nm )f mn (t). (34) 

Consequently, the s-population in cloud A in the point 
particle limit is given by 



neAmeB 



(35) 



Let us first consider the simplest case Na = 1 and Nb 
1 to determine the probability to find the s-excitation at 



time t on the atom in cloud A. Eq. (34) can be solved to 
give non-dephasing Rabi oscillations 

Pk(t) = |(7r 12 |*(<))| 2 = |/ 12 (t)| 2 = cos 2 (^] , (36) 



where dk denotes the distance between the two atoms, 
in a particular stochastic realisation k of the atomic spa- 
tial distribution. Averaging this function over a large 
number G of realizations yields the expression for the 
' s-population for point-like, fixed atoms, 



1 G 



(37) 



fc=i 



The expression does not change for the general case 
Na > 1, N B > 1: From the block structure of the elec- 
tronic Hamiltonian, one can infer that the averaged result 
is exactly the same as in the case N A = 1 and Nb = 1, we 
simply have added to the average over positions an av- 
erage over blocks, which has exactly the same structure. 
In the limit G — > oo, the sum in the above expression 
becomes an integral, in which the distances dk follow the 
initial spatial distribution. This gives precisely the same 
expression as in Eq. (31). In Fig. [2j we demonstrate the 
agreement between the exact quantum mechanical result 
and the point particle limit, where the latter is also shown 
for a larger number of atoms (N = 20). 

In summary, due to the particular structure of H e \, 
the time evolution of electronic populations per cloud 
for arbitrary Na, Nb is identical to the average over 
an ensemble of systems with Na = Nb — 1. We ex- 
plicitly demonstrated this using a point particle model, 
and further argue in Appendix|B]that this behaviour will 
persist fully quantum mechanically. Note, however, that 
while the point particle approximation gives a correct 
description for the specific case studied here, it will in 
general give different results than a full quantum me- 
chanical treatment. As an example, consider a setup in 
which the two clouds are additionally confined in two 
strong harmonic traps with energy spacing AE 3> Vg/L 3 . 
In this case the kinetic energy term is by construction 
much larger than the dipole-dipole interaction. Then, in 
a quantum treatment, -P^(i) would undergo long coher- 
ent oscillations without dephasing, while an average over 
point particles would show fast dephasing as without the 
traps. 



IV. ATOMIC MOTION 

We now extend our calculations to timcscales on which 
the dipole-dipole-induced atomic motion becomes clearly 
visible. Although the motion of the atoms is restricted 
to one dimension along the separation of the clouds, 
the full many-body problem is still too complex to be 
treated fully quantum mechanically. Therefore, we use 
a quantum-classical hybrid approach, which has been 
successfully applied to similar systems [14H16) . namely 
Tully's fewest switching algorithm [46-48]. The method 
relies on a classical treatment of the atomic trajectories 
and a quantum mechanical treatment of the electronic de- 
grees of freedom. A detailed description can be found in 
the aforementioned references, and we only briefly sum- 
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FIG. 2: Color online. S-population in cloud A as a function 
of time, shown for Na = Nb = 1 fully quantum mechanically 
(P|(t), Eq. \Vi\ , solid black line), as an average over point 
particles (p(f)7Eq. (371, red dashed line, for Na = Nb = 1), 
and as an average over point particles for Na ~ Nb = 10 
(crosses). The following parameters were used: a = /im, 
L — 6 fim, Vo = 10 6 a.u., corresponding to Li atoms with the 
principal quantum number v ~ 30. .40, and M — 11000 a.u. 
(mass of 6 Li). 



marize the main aspects relevant for the present problem 
in Appendix |C| 

The results of the simulations are shown in Fig. [3] We 
observe that, although the initial Rydberg excitation is 
coherently shared by all atoms within one cloud, it is 
nevertheless just a single atom that is ejected from each 
cloud. The atoms are still in a coherent superposition 
where each of them has left its cloud with some proba- 
bility, but the ejected atoms always carry the excitation 
and an ensemble of ground-state atoms stays behind. In 
this sense the superatoms break apart when subject to 
dipole-dipole induced motion, since excited and ground- 
state atoms are now spatially separated. 

Another interesting aspect is that the resulting mo- 
tion is a superposition of repulsion and attraction. It 
is known P3] that for single Rydberg atoms interacting 
via dipole-dipole-forces, the direction of the force can be 
controlled via the initial electronic state; The same holds 
for our system, i.e. the choice of | ip™) in Eq. (12 1 de- 



termines the motional dynamics. For Vo > 0, the initial 
electronic states 



i 



V2N A N E 



(kmn) ± kmn)) (38) 



n£A, m£B 



correspond to purely attractive (+) and repulsive (-) mo- 
tion of the atom pair. The initial state chosen in Eq. ( [12"] ), 
describing two superatoms, is a superposition of these 
two states, and consequently leads to a superposition of 
two opposite motion directions. 

The presented results have interesting implications for 
quantum transport protocols based on single-atomic ag- 
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FIG. 3: Color online. Atomic motion induced by dipole- 
dipole interactions. The upper panel shows the atomic density 
as a function of time. The lower four panels show snapshots 
of the density at times t = 0, 2, 4, 5 /xs. The color indicates 
the probability to find any Rydberg excitation (|s) or |p}) at 
a given position. In this figure, we do not distinguish be- 
tween \s) and |p) excitations, since the dephasing oscillations 
of electronic population, as shown in Fig. [2] occur on the 
short timescale t ~ 0.2 fis and are therefore not shown here. 
We observe that the main bulk of the atoms remains at rest, 
with vanishing probability to find an excitation. A pair of 
single atoms is ejected, and the Rydberg excitation is entirely 
localized on those two atoms. The calculations are shown for 
Na = Nb = 3, for the same parameters as in Fig. [2] but with 
M = 12600 a.u., corresponding to the mass of 7 Li. 



gregates, such as adiabatic entanglement transport [151 
116) . Our findings show that even in an ensemble of 
atoms which share the Rydberg excitation, only one 
of them will be set in motion by dipole-dipole forces. 
This facilitates the experimental realization of transport 
schemes in [151 I16| . since there is apparently no need 
to isolate single atoms. As formally discussed in Ap- 
pendix |Bj our results can be extended to more than two 
clouds, due to the block structure of the full Hamiltonian, 
Eq. ([2]), and the corresponding time evolution operator 
U(t) = exp [— iHt\, A 'Gedankenexperiment' presenting 
the transport scheme of Refs. [El [16] both, with single 
atoms and with atoms clouds, is sketched in Fig. [4] If 
atom clouds are used one should consider that the trans- 
port efficiency may be affected by collisions of moving 
Rydberg atoms with resting ground state atoms for a 
large number of atoms per cloud. However, as estimated 
in Appendix ID] in an ultracold gas such processes can 
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FIG. 4: Color online. Sketch of the excitation transport 
scheme with atom clouds as sites (lower panel), compared 
to the one with single atom sites (upper panel) as demon- 
strated in Refs. [151 116] . Color codes are as in Fig. [I] All 
but one atom from each cloud remain at rest and do not par- 
ticipate in the dynamics, hence excitation and entanglement 
transport occur equivalently in both schemes. For the sake of 
clarity, the lower panel shows the dynamics for a well defined 
choice of initially excited atoms in each cloud, while the com- 
plete dynamics consists of a superposition where each atom 
in each cloud has the chance to be the only one from that 
cloud participating in the dynamics. 






x[|j,m 



FIG. 5: Color online. Comparison of the atom density during 
superatom break-up, obtained from a full quantum mechan- 



ical calculation using Eq. (171 (solid black line) and Tully's 
quantum-classical hybrid algorithm (red dashed line), (a) To- 
tal atom density as a function of time, (b)-(d): Snapshots at 
three different times, as indicated by white vertical lines in 
(a): ti = 1/xs (b), t 2 = 3 ^is (c), t 3 = 4.9 /is (d). The same 
parameters as in Fig.[2]were used. For the sake of clarity, the 
comparison is show n for an initially purely repulsive electronic 
state, cf. Eq. p8l. 



AVOIDING THE BREAK-UP: 
DRESSING 



RYDBERG 



In the previous section, we demonstrated that dipole- 
dipole interactions between superatoms do not lead to 
their collective motion. In this section, we will show that 
the situation is quite different for essentially the same 
setup but using Rydberg-dressed ground state atoms [38] . 

For each atom, we now take into account two hyperfinc 
ground states \g), \h) and two Rydberg states \s), \p). 
Using off-resonant dressing lasers fields, we selectively 
couple \g) to \s) and \h) to \p). The corresponding Rabi 
frequencies and detunings are denoted by fl s , Q p and A s , 
A p respectively, see Fig. [6] and Ref. [38] . For simplicity, 
we take here the Rabi frequencies and detunings to be 
the same for both transitions, 



(39) 



The dressing lasers lead to a small admixture of the ex- 
cited states to the ground states, giving rise to dressed 
states 



be neglected on the microsecond timescale. Before con- 
cluding this section, we also mention that the validity 
of the adopted quantum-classical hybrid approach was 
verified by comparing the time evolution of the spatial 
probability density with a full quantum calculation for 
the accessible case = = 2, as demonstrated in 
Fig. [5] 



\g) =tf(\g)+a\. 



■ Af(\h) + a\p)), (40) 



where a = f2/(2A) <C 1 is a dimensionless scaling param- 
eter and J\f a normalisation factor. As shown in Ref. |38j . 
the dressing lasers lead to an effective interaction Hamil- 
tonian 



H cS = V 12 \g) 1 \h) 2 (h\ 1 (g\, 



(41) 



one atom 



s> 



lg> 



Vsp 



As 



lis 



/s - pair 



Ap 



l>p 



IP> 



|h> 

h/p - pair 



1 only by swapping with a second atom 



FIG. 6: Color online. General scheme of the dressing mech- 
anism. Only one atom is shown. For a transition, a second 
one with the reverse process is needed. The effective^ inter- 
action between two atoms in dressed ground states \gh) and 
\hg) can be sketched as \gh) — > \sp) ~* \ps) — > \hg), where '— 
symbolizes a dressing laser-driven transition and '*•»' a state 
flip via dipole-dipole interactions. In total, the scheme allows 
for an interchange of ground states \gh) and \hg). 



between two atoms 1 and 2, with 



Vi- 



i - [v nm { Trim 



(42) 



Illustratively, one can think of the interaction (41 1 as 
arising from a three-step-process: In the first step, the 
ground states \g), \h) are excited to the Rydberg states 
I s ): \p) by the dressing laser. In the second step, the 
dipole-dipole interactions flip the Rydberg states | s ) and 
\p). In the final step, the Rydberg states are de-excited 
back to the respective ground states by the dressing field 
(see Fig. |6| . This involves four photoinduced transitions 
in total, hence the interaction strength scales as a 4 . As 
a«l, the interaction strength between dressed ground 
states is smaller than between pure excited states. How- 
ever, the lifetime of dressed Rydberg states is larger than 
the one of actual Rydberg atoms [55] , 

Next, we will analyze these dressing- mediated inter- 
actions between two blocked clouds. We can distin- 
guish two regimes in Eq. (42 1: For separations much 
smaller than the dipole-dipole blockade radius rdd ~ 
(Vq/2H\A\) 1 ^ 3 , the interactions are strongly suppressed. 
For separations much larger than r^d we can simplify 



V 12 



a 4 V r . 



For cloud widths a < r^d and separations 

„4t 



reside in the same cloud, and V12 = otherwise. 

In the dressed case, we also have to modify the block- 
ade radius for Van-der-Waals interactions, which pre- 
vent the simultaneous excitation of two s or two p states 
within one cloud, to r V dw ~ (C 6 /2h\ A|) 1/6 [37]. The en- 
ergy level shift of doubly Rydberg excited states, caused 
by Van-der-Waals interactions, renders the dressing field 



ineffective, if this shift exceeds the detuning. Further 
care has to be taken that this offset and the detuning do 
not counteract each other, i.e. we assume a red-detuned 
dressing laser if the interactions are repulsive and a blue- 
detuned dressing laser if the interactions are attractive. 

Altogether we are led to an electronic Hamiltonian in 
the dressed case 



H*(R) = Vnm \g)n\h}rn(h\n(~g\r 



(43) 



which is a many-body version of (41) where no dressed 
interactions exist between atoms in the same cloud. This 
picture holds for cx 2 Na.b ^$ 1, otherwise our block- 
ade assumption becomes inconsistent with the mean ex- 
cited atom numbers a Na,b arising from the perturba- 
tive treatment in [35] . Since Eq. ( 43 1 conserves the to- 



tal number of atoms in either dressed one-particle state, 
let us restrict the following discussion to the case of 
even N with equal numbers of atoms in \g) and \h), 
as well as equal numbers of atoms in the two clouds: 
Na = N B = N s = N~ h = N/2. For the simplest non- 
trivial case N = 4 our basis contains the elements {\gg : 
hh),\gh : gh),\gh : hg),\hg : gh),\hg : hg),\hh : gg}}, 
where the 2 entries before the colon label the one-particle- 
state of each atom in cloud A, while the 2 entries after 
the colon label the states of the atoms in cloud B. For 
general N, the number of basis states is given by the 
binomial coefficient (N N/2). 

In terms of the basis above, the matrix elements of the 



dressed electronic Hamiltonian (Eq. 43 ) are given by 



Hi 



( 





V23 


V24 


v 13 


V u 





\ 




^23 














V u 






V24 














v 13 






v 13 














V24 






V u 














v 23 




V 





V u 


v 13 


V24 


v 23 





/ 



(44) 



In contrast to the Hamiltonian of bare dipole-dipole in- 
teractions in table |lj the one for dressed interactions (44) 
does not de-couple into blocks. 

Just as in the case of "standard" dipole-dipole inter- 
actions, the electronic Hamiltonian fl^(R) depends on 
the positions of the atoms, hence we again expect the 
interactions to give rise to motional dynamics. The full 
dressed Hamiltonian has the same form as in Eq. §}, 



N 



i=l 



2M 



Hi< 



R) 



(45) 



To infer the motional dynamics resulting from ( 45 ) , we 



consider the adiabatic energy surfaces defined in Eq. ( 14 ), 
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i.e. the eigenvalues of H^(R): 



Ui(R) 
U 2 (R) 


= o, 
= o, 










= J(y, A - 


y AO / l 




V24) 2 


U A (R) 


= -\M 4 


- V23) 2 


+ (Vl3 


- V2A? 


U 5 (R) 


= \J(V U + 


v 23 f + 


(Vl3 + 




U 6 (R) 


= -\/(Vl4 


+ v 23 y 


+ (Vl 3 


+ v 2i y 



(46) 



In a regime where the atomic motion is adiabatic with 
respect to the energy scale defined by ( 46 ) , the gradient 
of the k — th eigenvalue with respect to the position r.; of 
the i — th atom determines the force and hence motion 
of this atom IT4HT61. 



F? = -8 ri U k (R). 



(47) 



Let us for clarity's sake consider the special case where 
the distances between all atom pairs with atoms in differ- 
ent clouds are the same, i.e. r 13 = r 14 = r 23 = r 2 4 — g. 
Then the gradient of the eigenvalues U5 and U§ has the 
structure 



( t \ 

-a 



(48) 



where £ = S^/^Vqo 4 / q a is the magnitude of the force. 
We see that the atoms in cloud A both experience the 
same force, as do both atoms in cloud B (with opposite 
sign). In other words, the dynamics on these two adi- 
abatic surfaces corresponds to repulsion and attraction 
of the whole cloud. In the general case where the inter- 
atomic distances are not exactly the same (hence neither 
are the matrix elements Vtj) the magnitudes of the forces 
will slightly differ, which does not qualitatively change 
the picture as long as the spatial extension of each cloud 
is much smaller than the distance between the clouds. 
It is straightforward to extend the special case consid- 
ered here to larger atom numbers, i.e. one can verify that 
there always exists an adiabatic surface that corresponds 
to motion of the whole cloud. This statement is cor- 
roborated in Fig. [7J where we show the time-dependence 
of the total atomic density for the case Na = Nb = 4 
calculated with Tally's algorithm, starting in a wholly re- 
pulsive electronic state. It is also demonstrated that the 
dynamics is essentially adiabatic, with a non-adiabatic 
population transfer to other surfaces of the order of 10~ 7 . 

It is quite remarkable that the observed motional 
dynamics resembles collective behaviour, although the 
atoms within one cloud are non-interacting. This is due 
to the fact that all atoms in one cloud are initially lo- 
cated very close to each other and experience (almost) 
the same force. 




FIG. 7: Color online. Motion of dressed ground state atoms 
in a fully repulsive exciton state, (a): Total atom density as 
a function of time. Contrary to a superatom, the entire cloud 
of Rydberg-dressed ground state atoms is set in motion, (b) : 
Non-adiabatic population transfer from the adiabatic surface 
on which the dynamics is initiated is of the order of 10 -7 . The 
dynamics is shown for 8 atoms in total, with Na = Nb = 4. 
The following parameters were used: M = 12800 a. u., Vq = 
1.08 ■ 10 6 a.u., a = 0.5 fim, L = 8 fim, = 4 MHz, a = 0.18. 



VI. SUMMARY AND CONCLUSIONS 

We studied the motion of two Rydberg-blockaded atom 
clouds, induced by their mutual resonant dipole-dipole 
interactions. On timescales where the kinetic energy of 
the atoms is still negligible, the dynamics in the electronic 
subspace exhibits superexponential dephasing. This be- 
haviour can be very well described and understood within 
a mixed quantum classical framework, treating the nuclei 
as classical point particles and the electronic subspace 
as a discrete quantum system. For timescales on which 
the atomic motion is clearly visible, we observe that an 
initial coherent single excitation on each of the clouds 
(superatomic state) is physically removed by the dipole- 
dipole interactions. A single atom pair leaves the clouds, 
with the excitation being entirely localized on this pair. 
This finding may facilitate an experimental realization 
of quantum transport protocols in Rydberg aggregates, 
such as proposed in |13l 115) since one-atom-sites may be 
replaced by micro traps containing several atoms. 

In the second part of the paper, we considered a mod- 
ified setup where genuine Rydberg excitations were re- 
placed by a weak Rydberg state admixture to ground 
states. In contrast to fragmenting superatoms, the laser- 
dressed ground state atoms remain grouped together, 
and the motion from dressed dipole-dipole interactions 
is seen for the whole cloud. This opens a new direction 
for the implementation of motional mesoscopic nonclassi- 
cal states. Consider for example a scenario in which two 
such clouds are prepared in a superposition of collective 
attraction and repulsion. If interference upon recombina- 
tion of the two quantum pathways could be observed, one 
would have demonstrated a quantum mechanical super- 
position which extends over several micrometers. Such 
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a scheme presents a realization of an optically resolvable 
Schrodinger cat [49] . 



Appendix B: Atomic motion from block-diagonal 
electronic Hamiltonians 



Appendix A: Herman-Kluk Propagator 



In this appendix we simplify the Herman-Kluk propa- 



gator, Eq. (19), to arrive at the expression in Eq. (24), 



and elaborate on the conditions under which this simpli- 
fication is valid. 

The semi-classical trajectories encapsulating the frozen 
gas approximation are given by 



q« = q , 



Pt = p 



dU(q') 
dq' 



t. 



(Al) 



where q',p' are the initial phase space coordinates at 
t' = 0. Physically, this means that the atoms are too 
heavy to move on the timescale of interest, but there 
is nonetheless a build-up of momentum, which affects 
the phase. Let us first consider the prefactor C in the 
propagator (Eq. (21)). From the trajectories, we find 
the derivatives dpT/op' = 1, dq t /dq' = 1, dq t /dp' = 0, 



and dp t / dq' 

becomes 



-d U {q')/dq' . The prefactor therefore 



C = 



a 2 a 2 CT ± (q / ) 



1/2 



1± 



6V a 2 t 
ihL 5 



1/2 



ih dq' 2 

(A2) 

where we used the adiabatic surfaces U ± (q') = ±Vb/q' 3 
and set |q'| w L after taking the derivative. On short 
timescales (cf. Eq. (22)), we can thus approximate the 
prefactor by unity. Furthermore, since the classical ki- 
netic energy is, for the parameters of our setup on the 
given timescale, much smaller than the classical potential 
energy, we can approximate the action in the Herman- 
Kluk propagator as 

S ± [q t ,p t ]^-U ± (q / )t. (A3) 

Finally, we insert the trajectories into the explicit form 
of the coherent states, used for the propagation, and 
carry out the momentum integration. This yields the 
following expression for the propagator: 



K HK (R",t;R',t = oy- 



d N q' [ X %(R' - q')Y 



x exp 



C/^q') 



dq' 



(R"-q': 



<5(R'-R"), 
(A4) 



where we used the property of the delta function 
f(x)5(x — a) = f(a)S(x — a). The expression in the ex- 
ponent is nothing but the first two terms of the Taylor- 
expansion of J7 ± (R") around R" = q', which, in our 
case, is sufficient. Hence we replace the exponential by 
cxp(— itU ± {T{!')/K). The remaining integration over q' 
gives unity, according to the normalization of Xn- We 
thus arrive at the expression from Eq. ( 24 ) . 



In Section |HIj we have seen that the electronic Hamil- 
tonian for a collection of Rydberg-blockaded atom clouds, 
Eq. (pi), can be cast into a block-structure, where in each 
block only a single atom per cloud participates in non- 
trivial electronic dynamics. In this appendix we formally 
show that this leads to rapid entanglement of motion and 
electronic state, since in each block also only one atom 
per cloud performs non-trivial motional dynamics. This 
holds for any number of clouds and atoms per cloud. Fur- 
ther we show that non-adiabatic transitions do not affect 
this picture, as they respect the block structure. Let 
us consider the example of three clouds with two atoms 
each. The generalisation to other numbers is straightfor- 
ward. 

Let the atoms be grouped into clouds A = {1, 2}, B = 
{3,4}, C = {5,6}. In that case the total Hamiltonian, 
including atomic motion and the electronic degrees of 
freedom, can be written as 



H 



E 



T n + diag[H u H 2 ,H 3 ,...]. 



(Bl) 



We abrreviated the kinetic energy term for each atom 
as T n = -h 2 V 2 rn /2M. The symbol diag[ifi, H 2 , H 3 , ...] 
denotes a block-diagonal matrix, with 3x3 blocks 
Hj(rk,ri,r m ). Crucially each block only depends on the 
coordinates of three atoms, with k € A, I G B, m e C, 
i.e. one from each cloud. The latin index j just num- 
bers the blocks. Consequently, the time-evolution opera- 
tor U(0,t) = exp [—iHt/h] also obtains a block-structure 
U = diag[U 1 ,U 2 ,U 3 ,...], with 

Uj = exp[-i( T n)t/K\ 

x exp [-i{T k + Ti + T m + H 3 {r k ,r u r m ))t/h]. (B2) 

We have used the operator hat on U here exclusively 
to avoid confusion with the potential energy surfaces U 
occuring elsewhere. 

In Eq. ( B2 ) , the second exponential acts only on atoms 
k, I, m, while the pre-factor describes free motion of the 
remaining atoms. Now consider a time evolution begin- 
ning in an eigenstate of the electronic Hamiltonian H c \ 
which will have support only in a single block jo of H . 
The form of Uj immediately tells us that three particles 
undergo dipole-dipole dynamics and the other ones per- 
form free motion. Starting in a non-eigenstate will lead 
to quick entanglement of motion and electronic state, as 
in each block three j different particles participate in 
non-trivial dynamics. The block structure of Uj also im- 
mediately confirms the absence of non-adiabatic transi- 
tions between exciton states in different blocks. For an 
alternative argument, note that V H c \ will have the same 
block-diagonal structure as H e i. From that we can see 
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that non-adiabatic coupling terms, denned in the next 



appendix [Cj in Eq. C3 vanish whenever the involved adi- 
abatic eigenstates have support in different blocks. 



Appendix C: Tully's algorithm 

We start from the total Hamiltonian as given in 
Eqs. ([2]) or (Bl|. Tully's algorithm is implemented as 



follows. One first finds the adiabatic eigenstates in the 
electronic subspace ^(R)) and the corresponding ener- 
gies U k (R) , and expands the wave function in this basis, 
cf. Eqs. (Tl4|) and (15). However, contrary to a full quan- 



tum mechanical approach, the expansion coefficients do 
not depend on the atomic positions, and here we denote 
them by c k (t) to underline the difference: 



2N A N B 

*(R,t)= Yl c k (t)\<p k (R)). 
fc=i 



(CI) 



The adiabatic expansion is inserted into the time- 
dependent Schrodinger equation, which leads to a set of 
equations for the time-dependent expansion coefficients: 



2N A N B 

ihc k = U k c h ; - ih R-dfejCj. 

3=1 



(C2) 



In the above equation, we introduced non-adiabatic cou- 
pling vectors 

d kj = < W |V RW ) = Uj{B) _ Uk{R) ■ (C3) 

Atomic motion is treated classically, i.e. the time evolu- 
tion of the spatial degrees of freedom is obtained from a 
classical equation of motion (CEOM) 



MR = -V R E/ fe (R) 



(C4) 



and averaged over many trajectories. The averaging is 
performed such that the initial conditions for the CEOM, 
namely the position R(f = 0), the velocity R(t = 0) and 
the adiabatic surface k(t — 0), resemble the Wigner dis- 
tribution of the initial state. Then each trajectory is 
propagated on a single adiabatic surface k and the pres- 
ence of non-adiabatic couplings d k j is accounted for by 
introducing the possibility for instantaneous stochastic 
switches to another adiabatic surface j, as described in 
detail in Ref. l47l. 



Appendix D: Rydberg atoms moving through 
background gas 



corroborate this claim, in what follows we estimate the 
impact of the processes that could potentially upset a 
dynamics such as sketched in Fig. [4] Let us consider Al- 
kali atoms in an ultra-cold gas. They can be routinely 
excited to Rydberg states, and in general collisions with 
surrounding atoms do not necessarily lead to a loss of the 
excited state 50- 52J . However, in the setup considered in 
Section[TV] some Rydberg atoms are additionally acceler- 
ated within the gas due to the dipole-dipole interactions, 
with typical velocities of the order of v m = 4m/s for the 
example of 7 Li. This is about two orders of magnitude 
larger than the thermal atomic velocity in an ultra-cold 
7 Li gas at T = \\iK of vt — 4.8cm/s. However, as 
we estimate below, this does not lead to new regimes of 
atomic collisions. We consider the impact of the following 
processes: (i) Inelastic collisions of Rydberg and ground 
state atoms, (ii) elastic collisions of Rydberg and ground 
state atoms, and (iii) collisional ionization of Rydberg 
atoms. 

(i) We can estimate the rate T q of Rydberg state 
quenching collisions with background atoms, using T q — 
vn<j q {v), where v is the velocity of the Rydberg atom 
with principal quantum number v, n is the density of 
background atoms and <J q (y) is the cross section for the 
process. From [53] we have, e.g. for Rubidium atoms, 
cr 9 (40) = 2 x lQ~ u cm 2 . At v = 4m/s and with den- 
sity n = lxlO 18 m the inverse quenching rate is 
t = 1/Tq = 125 fis, which is still larger than the Ryd- 
berg state lifetime r = 46/xs. Note that for a scheme as 
depicted in Fig.[3]the effective background density would 
be an average over the intra- and intercloud regions and 
thus much smaller than the cloud peak density. 

(ii) For elastic collisions, Ref. [S3] gives a cross section 
of the order of a c \ « 10 4 a.u. « 0.25 x 10~ 12 cm 2 , which is 
about a factor of 80 smaller than the one given above for 
inelastic collisions. Elastic collisions are therefore even 
less important for the considered range of parameters. 

(iii) In order to estimate the order of magnitude of the 
rate for collisional ionisation, we adopt an approximative 
formula derived by Lebedev |55j . For v 3> I, it reads, in 
atomic units, 



r ion = exp 

7T/Z 



1 



(Dl) 



where a c \ is the cross section for elastic scattering of a 
Rydberg electron on a ground state atom, /1 the reduced 
mass of the colliding system, T the temperature and 
^eff = Si being the quantum defect. In the ultracold 

regime, the low temperature of ~ l^K w 0.3 x 10 -11 a.u. 
leads to a vanishingly small exponent, making the colli- 
sional ionization absolutely negligible. 



In Section |IV| we argue that excitation and momen- 
tum transport in Rydberg aggregates, as described e.g. 
in Refs. [THUS], should, in principle, be equally possi- 
ble if single atom sites are replaced by atom clouds. To 
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